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Abstract 



The dynamics of a single fluid bilayer membrane in an external hydrody- 
namic flow field is considered. The deterministic equation of motion for the 
configuration is derived taking into account both viscous dissipation in the 
surrounding liquid and local incompressibility of the membrane. For quasi- 
spherical vesicles in shear flow, thermal fluctuations can be incorporated in 
a Langevin-type equation of motion for the deformation amplitudes. The 
solution to this equation shows an overdamped oscillatory approach to a sta- 
tionary tanktreading shape. Inclination angle and ellipiticity of the contour 
are determined as a function of excess area and shear rate. Comparisons to 
numerical results and experiments are discussed. 

PACS Numbers: 47.15Gf, 68.15+e, 82.70-y 
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I. INTRODUCTION 



The equilibrium behavior of fluid vesicles made up of a single lipid bilayer membrane has 
been studied in great detail and is by now well understood . These equilibrium phenomena 
comprise the occurrence of a multitude of vesicle shapes and their thermal fluctuations. If 
external parameters such as temperature or osmotic conditions are changed, such a shape 
may become unstable and settle into the next dynamically accessible minimum of bending 
energy. The dynamics of such a shape change, however, is a non-equilibrium process. In 
general, only very little is known systematically and quantitatively about this and other 
non-equilibrium phenomena of fluid vesicles. 

A rough classification of non-equilibrium phenomena should distinguish between two 
types. First, there is relaxation into a new equilibrium after a parameter change. This 
class comprises the above mentioned decay of a metastable shape such, e.g., as the budding 
process |2],|3| . The best studied case of this class is the spectacular pearling instability of 
cylindrical vesicles, which develops upon action of optical tweezers [|]]. 

The second class of nonequilibrium behavior refers to genuine non-equilibrium states 
induced by external fields such as hydrodynamic flow fields. The behavior of vesicles un- 
der such conditions has only recently been started to be explored theoretically and 
in experiment ||. We developed a numerical code which allows to follow the shape evolu- 
tion of vesicle in shear flow |J. As a main result, we found that the shape finally settles 
into a non-axisymmetric stationary ellipsoidal shape around which the membrane exhibits 
tanktreading motion. Quantitative predictions include the inclination angle and the tank- 
treading frequency as a function of reduced volume and shear rate. This work left aside two 
interesting aspects. First, the regime of very small shear rates was difficult to reach within 
this numerical approach because of long relaxation times. Even though the numerical data 
seemed to indicate that shear is a singular perturbation, no clear assessment of this regime 
was possible. Secondly, fluctuations, which would be particularly important at small shear 
rate, were not included into the numerical algorithm. 

The purpose of this paper is twofold. First, we present somewhat more explicity the un- 
derlying continuum theory for the dynamical evolution of an incompressible fluid membrane 
in an external flow field thereby generalizing previous work on dynamics of membranes in 
quiescent fluids M and within a Rouse model [TTJ. The resulting equations are strongly 
non- linear and must, in general, be solved numerically. Analytical progress, however, is pos- 
sible for so-called quasispherical vesicles Jl 1 {T?J for which deviations from a sphere remain 
small throughout the parameter space. Focussing on such vesicles, it becomes possible to 
include the effect of thermal fluctuations which is the second objective of this paper. A 
main effect of thermal fluctuations in fluid membranes is that they store or "hide" area in 



the suboptical range [|18j. The present work addresses the issue of how external flow pulls 
out hidden area from thermal fluctuations. Apart from its fundamental significance such a 
study is motivated also by the so far only published experiments on the behavior of vesicles 
in shear flow ||. In this experiment initially spherical vesicles elongate under shear thereby 
pulling out area from thermal fluctuations. The theory developed here should be used to 
analyze such experiments. 

For a somewhat broader perspective, it may be important to point out that there ex- 
ist a vast literature on the behavior of a "soft sphere" in shear flow ranging from liquid 
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droplets fl9||20|j to elastic capsules |2T 
|27l or synthetic microcapsules [28 



3[ or inert ellipsoids [p4H26fl modeling red-blood-cells 
Many of these systems show ellipsoidal deformations 
with revolving or tanktreading motion. The details of deformation, inclination angle and 
tanktreading frequency, however, depend crucially on both the constitutive relation for the 
elasticity of the "interface" of the soft sphere and the dissipative mechanisms involved. In 
our case, the membrane is determined by bending elasticity and local incompressibility. The 
main dissipation occurs in the surrounding fluid. 

How thermal fluctuations are affected by shear flow seems not to have been touched 
upon in these studies of soft spheres. For the topologically much simpler case of on average 
planar membranes in a dilute lamellar phase, the layer spacing is predicted to decrease with 
increasing shear due to the suppression of fluctuations p9jp0 . 



This paper consists of five sections, of which this is the first. Section II discusses the 
physical model underlying this work and derives the fundamental equation of motion for an 
incompressible membrane in external flow. In Section III, we specialize to quasi-spherical 
vesicles in linear flow. In Section IV, we discuss in detail shear flow for which we solve 
the equations explicitly and compare to available numerical and experimental work. A 
summarizing discussion is presented in Section V. 



II. FORMALISM 
A. Physical model 



The starting point for the dynamics of a fluid membrane is the bending energy |3T 



F K = ^JdA(2H-C ) 2 . (1) 

Here, k is the bending rigidity, H is the local mean curvature and Co a spontaneous curvature 
in the case of intrinsically asymmetric monolayers or an asymmetric liquid environment. We 
neglect Gaussian curvature energy which would enter only through boundary values, since 
we focus on the dynamics of a closed membrane. For simplicity, we also neglect the fact 
that for closed membranes a second energy term becomes relevant which takes into account 
that the bilayer consists of two tightly coupled monolayers which do not exchange molecules 
[32|. 

The bending energy is the driving force for dynamical changes. Dynamics in the micron 
world of vesicles is overdamped, i.e. inertial effects can safely be ignored as can easily be 
checked a posteriori by calculating the corresponding Reynolds number ||33|| . Dissipation 



takes place both in the surrounding liquid and in the membrane, in principle. For giant 
vesicles of micron size, the dominant dissipation is viscous dissipation in the embedding 
fluid [j3"4] ) j3"5|| . Therefore, a full treatment of the hydrodynamics of this fluid is mandatory for 
a faithful description of the dynamics of membranes. 

Dissipation in the membrane can be classified into three phenomena: Drag between 
the two monolayers ||36|| , shear viscocity within each layer and permeation through the 
membrane. Calculation of the relaxation spectra of bending fluctuations involving the first 
two mechanisms show that on scales of microns and larger, hydrodynamic dissipation is 
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dominant |37| . In the submicron range, friction between the layers becomes relevant. On 
even smaller scales of several tenths of nanometers, shear viscosity within each layer should 
be included. Finally, permeation through the membrane seems to be irrelevant on all length- 
scales with the possible exception of membranes in the vicinity of a substrate 



Based on this hierarchy of dissipative mechanisms, we will include only hydrodynamic 
dissipation and model the membrane as impermeable for liquid. As a consequence, the 
normal velocity of the fluid at the membrane pushes along the membrane and leads to a 
shape change. Tangential motion of the fluid along the membrane induces lipid flow within 
the membrane since we employ non-slip boundary conditions between fluid and membrane. 
One could allow for some slip with a phenomenological friction coefficient. In the absence of 
any evidence for such a phenomenon, however, we choose for simplicity the no-slip condition 
used so far. Finally, we require that the membrane remains locally incompressible. 



B. Elementary differential geometry of the membrane 

The instantaneous membrane configuration R(s 1 ,s 2 ) ) parametrized by internal coordi- 
nates (si, sa), is embedded in the three dimensional space. This space will be parametrized 
by Cartesian (x, y, z) or spherical (r, 6, <fi) coordinates where a = x,y, z. 

Summation over double indices is implied throughout the paper. There are two tangential 
vectors 

R t = <9iR(si, s 2 ) for i = s l5 s 2 , (2) 
from which one obtains the metric tensor 

9ij = Ri • Rj- (3) 
Its determinant, g = det^), yields the area element 

dA = y/g ds\ds2- (4) 
The normal vector, n(si, s 2 ), is given by 

Ri x R2 

Finally, the mean and the Gaussian curvature follow from the curvature tensor 

hij = (didjR) ■ n (6) 

as 

H = \h\ (7) 

and 

K = det{h l X (8) 
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where hj = g lk hkj and are the matrix elements of the matrix inverse of (gij)- Following 
the convention used in differential geometry, a sphere with the usual spherical coordinates 
(si = 8, s 2 = 0) has H < 0. A membrane configuration has bending energy F K given in ([I]). 

In order to ensure local incompressibility of the membrane, we will need a local Lagrange 
multiplier £(si, s 2 ) which we call local surface tension. It will be determined self-consistently 
below. The total "energy" thus becomes 

F = F K + / d Sl ds 2 ^g- E(si, 8 2 ). (9) 



C. Equation of motion 

Any non-equilibrium membrane configuration exerts a local three dimensional force den- 
sity 

f ( r ) = -J dsi ds 2 ^ (-^) S(v - R(s h s 2 )) (10) 

onto the surrounding fluid. 

The variational derivative entering the force density reads explicitly 

(--)- 



-2EH + k[(2H + C )(2H 2 -2K- C H) + 2AH]) n - ^R^S. (11) 



Here, 

A = (l/yffiditfty/gdj) = g'>»,» r (12) 
is the Laplace-Beltrami operator on the surface. The normal part of (|TTD is well-known from 



the stationarity condition of membrane configurations | 39| . The tangential part arises from 
inhomogeneities in the surface tension which will be needed to ensure local incompressibility 
of the induced flow. This approach differs from the one used in Ref. where a finite 
compressibility was introduced for renormalization purposes. 
The surrounding liquid is incompressible 

Vv = (13) 

and obeys the Stokes equations 

Vp-7jV 2 v = f(r). (14) 
A special solution of the inhomogeneous Stokes equation for the velocity field reads 

v ind (r) = f dr'e>(r,r , )f(r'), (15) 



where the Oseen tensor 0(r, r') has Cartesian matrix elements |40 
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o i / 1 I — / 

8^77 |r — r | [ |r — r 



O a0 {r,r ) = ^—j- — — 5 af3 + - — . (16) 



Thus, the hydrodynamics mediates a long-range interaction (~ l/|r — r'|) through the ve- 
locity field. To this velocity field v ind (r) induced by the presence of the membrane, we must 
add the externally applied flow field v°°(r) to obtain the total velocity field 

v(r) = v°°(r) +v ind (r). (17) 

Such a simple superposition becomes possible since the Stokes equation is linear due to the 
absence of the convective term. The value of the total velocity field at the position of the 
membrane yields the equation of motion for the membrane configuration 

d t R( Sl , s 2 ) = v°°(R(s 1 , s 2 )) + v ind (R( Sl , s 2 )), (18) 

since we employ no-slip boundary conditions at the membrane. This deterministic equation 
of motion includes both normal motion that signifies a shape change of the membrane and 
tangential motion that corresponds to lipid flow within the membrane. The so far unknown 
local tension S(si,S2) is determined by requiring that the membrane flow induced by this 
equation obeys local incompressibility. Thus, we must demand 

d t y/g = 0, (19) 

which implies explicitly 

gVRidjidtR) = 0. (20) 

Upon insertion of the equation of motion (18) with (|T0|) and (|T5|) , this condition becomes a 
partial differential equation for the unknown tension E(si, Sg). 

Equations fllS]) and ( pO"D yield a deterministic evolution equation for a membrane con- 
figuration under the action of bending energy and hydrodynamics. For any given initial 
membrane configuration, the solution to these equations will run into the next dynamically 
accessible local minimum of bending energy. In general, these equations must be solved 
numerically. For vesicles in shear flow, this has been achieved recently |J. 

III. QUASI-SPHERICAL VESICLE IN FLOW 



In this section, we apply the general formalism to quasi-spherical vesicles fIT|-|T7|j for 
which analytical progress becomes possible. 

A. Expansion around the sphere 

A quasispherical vesicle is characterized by its volume 

V EE f R \ (21) 
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which defines R, and its fixed area 



A = (4tt + A)R 2 , (22) 

which defines the (dimensionless) excess area A. 

The instantaneous vesicle shape R(6 I , 0) = R(6, 0)e r can be parametrized by spherical 
harmonics 

(lmax I \ 

X +E E <WWM) , (23) 

where \m\ < I and ui^ m = {—l) m u* lm . The upper cutoff l max is of order R/d where d is 
the membrane thickness. For a vesicle with R ~ 10 — 50fim, l max is of order 10 4 . Since 



spontaneous curvature is irrelevant for the shape behaviour of quasi-spherical vesicles [17 



we set for the rest of the paper Cq = 0. Expanding the geometrical quantities as well as the 
bending energy around a sphere, one has [|i^ , |i~3l . |3T| 

1 lmax I 

Fn = 8tt + - E E M 2 ( l + 2 )( l + W - 1) + 0«J, (24) 

/ i>0 m=-Z 



A = R z 



4 T + tIJ + it ^j^i 2 ( i+/ ( /+i )/ 2 ) + 



(25) 



and 



4tt A mq.o 
3 \ V4¥ 



3 L 



+ E E k/, m | 2 + Of 



(26) 



1>1 m =-l 

The volume constraint ( pT|) fixes the amplitude tio,o as a function of the other amplitudes 



lmax I 

u o,o = - E E KmlVv 7 !^, 

Z>1 m =-/ 



(27) 



where we truncate from now on the cubic terms. If this value is inserted into (|25|) , the area 
constraint (E2|) becomes 



; ' 2 (/ + 2)(/-l) 



E E k-. 

1>1 m =-l 



A 



(28) 



Since the (I = l)-modes correspond to translations, which have to be omitted, from now on 
all sums start at / = 2. We will abbreviate 



lmax I 

E= E E ■ 

Lm l>2 m=—l 



(29) 
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We now add the global area constraint (p8|) with a Lagrangean multiplier 

£ = na/R 2 



(30) 



to the quadratic part of the curvature energy. This leads to a quadratic expression for the 
total "energy" (|) 

F = ^5>,K m | 2 , (31) 

l,m 

with 

E t = (I + 2)(l - + 1) + a}. (32) 

The instantaneous state of a vesicle must be characterized not only by the set of its 
deformations u^ m but also by its instantaneous surface tension 

E(M) = ^(* + E<WWM))- (33) 

The homogeneous value cr = o"o,o, which will be called effective tension, has already been 
included into the energy (5^). 



B. Velocity field 

The notion of a quasi-spherical vesicle implies that the deviations u(9, </>) from the spher- 
ical shape are small. Then stresses caused by bending moments and an inhomogeneous 
surface tension can be assumed to act on the sphere rather than on the deformed vesicles's 
surface. Likewise, all velocity fields will be evaluated and matched at the sphere. We retain 
this procedure for arbitrary flow strength even though it is strictly valid only for small ex- 
ternal fields. As will be discussed in Sect. IV. D below, while not modifying scaling laws, 
this approximation causes an error in numerical prefactors at large shear rates of at most 
30%. 

The instantaneous stress distribution caused by both the deformation and the inhomo- 
geneous surface tension induces an irrotational velocity field v md (r) which we determine in 
the Appendix adapting the classical Lamb solution f33f . For this purpose, it is convenient 
to define another three-dimensional velocity field 

V ind (r) = v md {R,6,<j)) (34) 

which extends the velocity field on the sphere corresponding to the vesicle formally to all 
space. V md (r) is r-independent. This velocity field can then be characterized by its normal 
component 

X md _ £ xi^y^e, <P) = v ind (i?, e, <p)e r (35) 

l,rn 

and the negative of its divergence 
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^ ind = E Yi£yi,mW> <t>) = ~ R VV M , (36) 

l,m 

where y^ m {9,4>) are the usual spherical harmonics. Although the full three-dimensional 
velocity field v ind (r) is divergenceless, the derived quantity V md (r), in general, is not. 

Normal and tangential stress balances on the vesicle surface determine v md (r) as a func- 
tion of the displacement u(9, 4>) and the surface tension a(9, <p). As shown in the Appendix, 
we can determine the linear relationship between the expansion coefficients of the induced 
velocity on the sphere {Xf^, Y^} and the sources of this field {iiz, m , cx/ im }. Replacing a^ m 
by Xf^ and Y^, one finally gets 

4m = ~(k/vR 2 ) riEiU l>m - BiY™, (37) 

where 

Tl = 4/3 + 6/2-1 (38) 

and 

2/ + 1 

It is convenient to characterize the external flow field v°°(r) in a similar manner. We 
therefore define 

V°°(r)=v 00 OR,M), (40) 

its normal component 

X°° = E X? m yi,m(0, 4>) = V°°(/2, 9, <P)e r (41) 

l,m 

and its divergence 

^°° = E 4) = - W°°. (42) 

l,m 

The equation of motion (|T8| ) for the quasi-spherical vesicle in this external flow thus 
becomes 

d t R(9, <P) = V°°(fl, </>) + V ind (fl, 0). (43) 
The local incompressiblility condition ( f20| ) becomes after little algebra 

V(V°° + V ind ) = 0, (44) 

or 

yind = _yco ^ 
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Physically, the incompressibility condition fixes the local surface tension which we have 
already eliminated in favor of Y md . 

We will now set up the equation of motion for the normal deviation, u(8,<f)), from the 
perfect sphere, where 9 and are fixed coordinates in the lab frame. It is important to realize 
that (|43| ) leads to a tangential motion of the material point labeled by (#, 0). Therefore it 
is not sufficent to project ( |43D onto e r in order to get the normal shape change in the lab 
frame. We also have to include an advection term as for planar membranes [2^ , |5D ]. The 
equation of motion thus becomes 

d t u{6, <p) = -(V°° + V ind )Vu(9, cj)) + (X°° + X' md )/R. (46) 

Note that the advection term can be written in terms of three-dimensional vectors because 
u carries no r dependence. In general, the advection term will couple the different modes. 
Further progress becomes possible for the quite wide class of linear external flow. 



C. Langevin equation for linear flow 

Linear external flow is characterized by the form 

v°°(r) = Qr (47) 

where Q is a traceless r-independent matrix. It can be decomposed according to 

v°°(r) =Q s r + n x r, (48) 

where ft denotes magnitude and direction of the rotational part of the flow and Q s is a 
symmetric traceless matrix describing the strain or irrotational component of the flow. The 
strain component of the external flow is compensated by the corresponding induced flow 
because of the incompressibility constraint. The advection term thus involves only the 
rotational component of the flow and reads 

-(fix r)Vu(0, 0) = -in(L/7t)u(0, 0), (49) 

where L = r x (— «7iV) is the angular momentum operator. Choosing co-ordinates such that 
O = Qe z , we finally obtain from (|46|) the equation of motion 

d t u l>m = -zfim u l>m - (s/^jr,^ Ul , m + (Xr m + BiY^/R + (i,m (50) 



after expanding in spherical harmonics and using (|37j) and ( flop 

In order to determine the correlations of the thermal noise <^,m we apply this equation of 
motion to equilibrium [pTTfTl3|] . Then, there is no external flow, Q = Xf^ = = 0. In the 
long time limit, the dynamical equal time correlations calculated from fl50|) will reproduce 
the static correlations only if we choose 

< Ci,m(t)(i>,m>(t') >= 2T(i? 3 /r ? )r i (-l) m 5 M '^,- m '5(t - t'). (51) 
Here, T is the temperature and Boltzmann's constant is set to 1 throughout. 
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We keep the noise correlations (51]) in the presence of external flow. While this is certainly 
correct for small flow it is not clear whether and how strong flow modifies these correlations. 
Note in passing that we have refrained from adding noise to the general equation of motion 
(UP, since the appropriate correlations for these forces are not clear for the same reason. 
Naively, one would expect them to exhibit long-range spatial correlations given by the Oseen 
tensor for small flow. Moreover, even in equilibrium, there are subtleties associated both 
with the incompressibility constraint and measure factors |J. 

The stationary value u^ m of the amplitude Ui >m (t) follow from (0) as 

(riR 2 \ ( X%, + B,Y& \ 

where 

fi = Qt}R 3 /k. (53) 

For linear flow, ui y7n ^ only for I = 2. 

The deviations from this stationary value, 

ei, m (t) = Ui tm (t) - Ui, m , (54) 

obey the equation of motion 

d t % m = -iSlm % m - (k/i]R 3 )TiEi e^ m + 0, m . (55) 
This simple equation is easily solved as 

e hm (t) = exp[(-( K /r/ J R 3 )r i ^ - iftm)t}( f drexp[((«/^ii! 3 )r^ + ^m)r]O lW (r) + ej, m (0)). 

J o 

(56) 

Starting from an initial value ej jOT (0), the approach to the stationary value ( p2|) thus happens 
via damped oscillations. 

Using the noise correlations ([51]), the dynamical correlation function in the long time 
limit become 

lim < e, >TO (t)e,,_ m (t + At) >= (-l) m exp[(-(/ t /r ?J R 3 )r^ i + imQ)At}^- (57) 
t-*oo nh/[ 

with the stationary correlations 

< e lim {t)e^ m (t) >= (-l) m ^r. (58) 
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D. Area constraint and effective tension 



As a last step, we have to eliminate the yet unknown Lagrange multiplier or effective 
tension a in favor of physical quantities using the area constraint fl2~8|) . Excess area is stored 
both in the systematic stationary amplitudes ui >m as well as in fluctuations e^ m . For the 
systematic part, we get 

A = £ — 2 — Kml = 2 1— J I nwrsw J • (59) 

The fluctuation contribution A/ m of the mode e/ jm to the excess area is 

Ai , m s (i±Mzi),^|. = £±^>|, (60) 

The total area constraint becomes 

A = 2^V V ( W^ + ^&H T (Z + 2)(Z-1) 

V « / n^-AnmW + tern*) 2k fa E l {a) ' 1 ) 

where we made the cx-dependence of E\ explicit. This master equation determines the 
unknown effective tension a as a function of the physical parameters T/k and those charac- 
terizing the external flow. In general, it must be solved numerically. 



IV. SHEAR FLOW 

In this section, we specialize the general results of the previous section to the experimen- 
tally important case of simple shear flow with shear rate 7. 

A. Flow parameters 

We choose the coordinate system such that the externally imposed simple shear flow 
reads 

v°°(r) = ^ye x = i[{y/2)e x + (x/2)e y )] - (7/2)6, x r. (62) 
We have separated the rotational component and identify 

tt = -7/2 and tt = -x/2, (63) 

where 

X = T7#7« (64) 
is the dimensional shear rate. The elongational first part of the shear flow can be written as 
j[(y/2)e x + (x/2)e y )] = 7r(sin 2 #sin0cos(/>e r + 
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(1/2) sin 9 (cos 2 <p — sin 2 0)e^ + sin 9 cos 9 sin cos </>eg) . (65) 

One thus has 

X°° = Y^ = ^sin 2 0sin0cos0 (66) 
and identifies the expansion coefficients in spherical harmonics as 

X™±2 = = T^(27t/15) 1 / 2 (67) 
and any other Xf^ = = 0. Note as an aside that X°° = Y°° holds for any linear flow. 

B. Effective tension 

We must eliminate the effective tension a in favor of the area constraint. The external 
shear flow implies the non-zero stationary amplitudes fl52| ) 



TO(12/11)(2tt/15) 1 / 2 
U2,±2 = =r= : ' ( 68 ) 

r 2 £ 2 =F IX 

where we used B 2 = 1/11. According to ( |59|) the systematic part of excess area stored in 
the modes W2,±2 contributes 

A = 2(\u 2 , 2 \ 2 + \u 2 ,^\ 2 ) = a 2 f (69) 

1 2-^2 + X 

where 

a 2 = 4(12/11) 2 2tt/15 ~ 1.994. (70) 

The area stored in the fluctuations becomes 

T ^ (Z + 2)(Z — 1) _ T ^ 21 + 1 T 5 T ^ 21 + 1 

2k f£ Ei ~ 2K^ 2 l{l + l)+a ~ 2k6 + cj + 2k f^ 3 l{l + 1) + a' ^ } 

Replacing the last sum by an integral, we obtain the total area constraint in the form 

° 2 rpfo + £ (^) + (72) 

This equation can easily be solved numerically for a = cr(x,T/K,l max ). However, it is more 
instructive to discuss limiting cases and the general behavior analytically. 

For vanishing shear rate, i.e., in equilibrium, the area constraint implies the three regimes 
0): 

(i) Tense regime: For A << T/2k, one obtains from (f72j) 

° W 2^A Z — • (73) 
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In this regime, all N modes share the available excess area. 

(ii) Entropic regime: For T/2k << A << (T/k) lnl max , the tension depends exponen- 
tially on the excess area |T3 



o - t ax e- 2KA/T . (74) 
(hi) Prolate regime: For (T/k) ln/^ ax << A <, 1, most of the excess area is stored in the 



(I = 2)-modes. The tension approaches the limiting value —6 flT3l , |T7 



^- 6 + 5^ < 75 > 

A vesicles that is relaxed with respect to its volume constraint also belongs to this regime 
with cr = 0. 

For large shear rate, i.e., formally x ~ ¥ °°> the whole excess area is stored in the w 2j ± 2 
deformation. From 

A = 2(|« 2i2 | 2 + K_ 2 | 2 ) = a 2 f - = A, (76) 



one obtains 



and hence, with E% ~ 4cr and r 2 = 6/55, the limiting behavior 

a w (55/24)a 2 /2 x/A 1/2 ~ 3.24X/A 1 / 2 . (78) 

Thus, for large shear rate, the tension increases linearly with shear rate with a prefactor 
that depends on the excess area. 

The cross-over between the equilibrium tension and this shear rate dominated saturation 
regime happens at a crossover \ c . The scaling behavior of Xc can be obtained by defining 
Xc as the value for which half of the excess area is stored in the systematic contribution A. 
This leads to the expression 

Xc ~ <i l 2 max e- 2KA/T \ entropic, (79) 
T/k A, prolate. 



C. Shape parameters 

Knowing the value of the effective tension as a function of the excess area, we can de- 
termine the characteristic parameters of the shear rate dependent shape of a quasi-spherical 
vesicle. 

The deviation of the mean contour from a circle in the plane z — 0, i.e., 9 = tc/2, becomes 
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«(tt/2,0) = 2^e{ M2i2 y 2 , 2 (vr/2,0)} = cos 2 (0 - O ), (80) 



with the mean incination angle 



- arctan — - — -. (81) 
2 X K } 



Thus, the mean incliniation angle is 0o — tt/4 for small shear flow. For large shear rate, the 
limiting behavior 

1 /a.-AX 1 / 2 , AV2 

0o « - arctan ^ — j « vr/4 - (82) 

follows from (|77|). Thus, the more excess area is available, the smaller the inclination angle. 
The cross-over between the two limiting cases happens at the crossover shear rate \ c given 
in (0). 

Experimentally, the ellipticity of such a contour is often measured by the deformation 
parameter 

n - L ~ B - 6 x 
= TTb ~ TiplEZ + x 2 ) 1 ' 2 ' 

where 2L and IB are the long and short axes, resp., of the contour. For small shear rate \i 
the equilibrium scaling of the tension ( 73f) - (|75|) in the three regimes implies the following 



linear behavior of the deformation parameter D on ^ as 

2 k A 



D ~ 6 X _. J? X 5 ; 

~ TTr 2 E 2 ~ 11 (24/55) (6 + 0") ~ 4 X ' ^ 



T/2 > tense > 

J '"max 



e 2KA/T /C entropic, (84) 
[ (2kA/5T), prolate. 



With increasing shear rate, deviations from this linear behavior become important at 
the crossover \c (FID- For X >> Xc> the deformation saturates at the value 



D « v /15A/32tt, (85) 
obtained by combining fl83|) with ( p5|) . 



D. Comparison to numerical work 

For large shear rate, fluctuations become irrelevant. In this regime, we can make contact 
with previous numerical work where the continuum equation of motions ([11]) and fl20"|) were 
solved on a triangulated mesh ||. The two quantities that can be compared within the 
two approaches are the dependence of the inclination angle and the effective tension on the 
excess area. In the previous work, reduced volume 
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v = , A ,owL M/o = (! + A/4vr)" 3/2 (86) 
(47r/3)(A/47r) 3 / 2 v 7 y v ' 

was used instead of the excess area. This definition implies 

A = 4tt{v~ 2/3 - 1). (87) 
Hence the present theory predicts for large shear rate and 1 — v « 1 the limiting behavior 

7r / 2n(l — v) \ ix 



r^~> I ; ' 1 r-^j 



" 4 V 3a 2 ; 4" L025 ( 1 - " )in m 

This square-root scaling is in quantitative agreement with the previous numerical solution 
0. The numerical prefactor within the present approach is about 10% larger than the one 
extracted from Fig. 2 of Ref. ||. Likewise, the scaling of the effective tension (|78|) agrees 
with the numerical data |4T| with a difference in prefactor of about 30%. 

The origin of the numerical difference in pref actors presumably lies in neglecting higher 
order couplings between deformation and external flow in the equation of motion fl50|). These 
effects would lead to non-diagonal and non-linear terms in the right hand side of ( p0|) which 
would spoil the solvability. It is therefore gratifying to know that at least for the inclination 
angle of the tanktreading state, the present treatment is quite reliable for v > 0.8 at any 
shear rate. 



E. Comparison to experiment 

There seems to be only one published experiment on quasi-spherical vesicles in shear flow 
||. In this work, two different types of behavior have been reported, (i) Initially spherical 
vesicles rotate at low shear rate and deform at higher shear rates into a tanktreading state 
with an inclination angle close to tt/A. (ii) Initially non-spherical vesicles undergo a periodic 
flipping motion at small shear rate. This flipping motion shows oscillations at larger shear 
rate. Only for sufficiently high shear rate, stationary behavior with tanktreading is observed. 

The first class of behavior fits well into the theoretical description developed here. The 
experimentally found D(j) curve shows linear behavior at small 7 and seems to saturate at 
larger shear rate, de Haas et al || analyzed this curve using a theoretical approach which 
may look superficially similar to the one developed here. For future reference, we discuss 
briefly the differences between these two approaches. De Haas et al use 

. 4ED 64™ - 

T= 5iSr P 15T D ' (89) 

in our notation. The corresponding expression within our approach derived from (^) reads 

4DS , Qk , ( fllD\ 2 Y 1 ' 2 , . 

+ )) (90) 

in dimensionalized units. The k term in the first parenthesis arises from keeping bending 
energy in the force balance . De Haas et al deliberately ignore bending as being irrelevant 
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for large tension which is true in the tense and upper part of the entropic regime. Then 
both expression yield the same linear behavior at small D. The non-linear regime, however, 
is markedly different with no exponential term in our theory. Note also that, within our 
approach, E still carries an implicit 7 dependence. 

The second type of reported behavior cannot so easily be reconciled with the present 
theory which predicts a stationary tanktreading shape for all values of shear rates and 
excess area. However, the theory allows for a transient damped oscillatory approach to this 
stationary state. In the prolate regime, the relaxation time 

r 2 = t]R 3 /kT 2 E 2 (91) 

can become large for large enough excess area because then formally E 2 — > as A — > 00. 
For large enough shear rate, the product of oscillation frequency 7 with the relaxation time 
r 2 can become large. In this case, many oscillations can occur before the stationary state is 
reached. Thus, after turning on shear flow or changing the shear rate one expects transient 
flipping motion and oscillatory behavior. Still, the present theory predicts finally a settling 
into a stationary tanktreading state. 

V. DISCUSSION 

For a fluid membrane in an external hydrodynamic field, we derived a deterministic 
equation of motion taking into account both viscous dissipation in the surrounding liquid 
and local incompressibility of the membrane. For quasi-spherical vesicles in linear external 
flow, this equation decouples in an expansion in spherical harmonics. The area constraint 
is implemented by an effective tension that depends on shear rate. Adding noise with 
correlations adapted from equilibrium, we obtain a Langevin-type equation. Its solution 
for shear flow yields transient oscillations which settle into a stationary tanktreading state. 
The inclination angle and the deformation parameter are both calculated as a function of 
shear rate and excess area. For small shear rate, the inclination angle becomes ir/4. With 
increasing shear rate it decreases the stronger the larger the excess area is. The deformation 
parameter increases linearly with shear rate with a slope depending on excess area. 

This work is complementary to our previous study where we solved the deterministic 
equations numerically for arbitrary excess area ignoring fluctuations ||. Taken together, 
the two approaches yield a fairly complete picture for the behavior of fluid vesicles in all 
phase space. There remains only one explored region. For an excess area too large for 
the quasi-spherical approximation to apply and a shear rate too small for the numerics to 
handle, neither approach is well-suited. In this regime, one would naively expect that small 
shear leads to small distortions around the fluctuating equilibrium shape. However, this 
remains to be checked by an explicit calculation. 

A relatively straightforward extension of the present theory would be the inclusion of fur- 
ther dissipative mechanisms such as inter-monolayer drag and shear within each monolayer. 
While both types will certainly affect the transient behavior, it is less likely that monolayer 
shear viscosity changes the stationary results for quasi-spherical vesicles since the rotational 
component of the flow does not lead to shear within the layer. Likewise, one could easily 
allow for different viscosities inside and outside the vesicle. 
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The analytical progress possible within the quasi-spherical approximation comes with two 
intrinsic weaknesses which cannot easily be circumvented. First, in the prolate regime, this 
approximation gives only a rough representation of the stationary ellipsoidal mode since it 
does not produce a non-zero -u 2 ,m without shear. Second, this approximation neglects higher 
order couplings between deformation and flow field which can become relevant at larger 
shear rates. However, the relatively good agreement between the present theory and the 
previous numerical work shows that this effect is certainly not dominant for v >, 0.8. 

The present study should encourage further experimental work. By measuring the in- 
clination angle and the deformation parameter of a single vesicle over the whole range of 
shear rates and fitting those against theory using the full equation flT2|), one should be able 
to extract fairly precise data on excess area and bending rigidity especially when combined 
with a variation of temperature |E3 . 

Concerning the experimentally reported oscillatory or flipping motion, it seems to early 
for a final assessment. While it could be that the above mentioned effects mask such a 
motion within the theory, there are two more indications for doubting that flipping motion 
is genuine for non-spherical vesicles. First, our numerical work || which did not suffer from 
the spherical approximation, did not show flipping motion for at least x ~ 0.1. Secondly, 
contrary to the reference given by de Haas et al ||, previous work on tanktreading of red- 
blood-cells [fl4] did not predict flipping motion if the viscosities inside and outside are the 
same. Only if the inner viscosity is significantly larger that the outer, a transition to flipping 
motion occurs. 
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APPENDIX A: ADAPTION OF LAMB'S SOLUTION 



Lamb's solution tells us the velocity field v(r) and pressure field p(r) in all space if the 
velocity field is specified on a sphere with radius R. We adapt the presentation of this 
classical problem given in Ref. |33| whose notation we follow. Rather than characterizing 
v(i?, 6, 4>) by its three components, Lamb's solution uses the three quantities 



XI 



= V(r)e r 



(Al) 



Y(9,(f>) = -R VV(r), (A2) 

and 

Z{9,<j>) = R e r (V x V(r)), (A3) 

derived from the field V(r) = v(R,8,<p). Note that V(r) is independent of r. Since the 
rotational component Z cannot be excited by either bending moments or inhomogeneities 
in the surface tension, we can ignore it in the following. 
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Given these boundary values of the velocity, the total velocity field for r < R becomes 

v-w = g (v*r(r) + ¥ ;;; ti) ^) - „ (f+1) ' (a+3) ^w) ^ 

where 

$r(r)= £ ^ >m ( M )( r /fl)i (A5) 

m=— I 

and 

Pi n (r)= E pjpWMXr/i?)'. (A6) 

rra=— £ 

The expansion coefficients are determined by the boundary conditions as 

^(') = §[(i + l)^l,m-nm] (A7) 

and 

Pfn» = ! ^^[^^ - ( / - ( A§ ) 

For r > R, the velocity field can be obtained formally by replacing I by — / — 1 in these 
expressions. Explicitly, it reads 

v~(r) = £ fv«f-(r) + ^^-^'Vrf-W - ^ ij rtf'w) (AS) 



Z>1 



where 



$r(r)= E S^ l>m (M)(r/i2) ( - ,+1) (A10) 



and 



/ 

m=— I 



Pi U \r) = E OWM)(r/i?) H+1) (All) 



with 



and 



p ^ (r) = Trri)§ [F/ ' m + (/ + 2)x '' m] - (A13) 
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This velocity field leads to a stress vector acting across the surface of a sphere of radius 
r = R as 

T R ee T R e r + T t ee (-e rP + ^ - X) + ^V(rv)) ^ . (A14) 

The normal component of this stress vector at r = R reads for the inner solution 

Tr=J2 {^/R?){1 - l)Z$f m + (Zof 1 - 6r)pf m ) M >m (*> 0), (A15) 



l.m 



where 



1(1 + 2) 



(/ + l)(2/ + 3) 



(A16) 



and 



b™ = 212 + 41 + 3 (A17) 
bl ~ (/ + l)(2/ + 3)- (A17) 

Likewise, the normal component of the stress vector at r — R reads for the outer solution 

T° R ut = e (Pfa/^x* + m + - [(/ + i)«r + &r>si) yuo, </>) ( Ai «) 

where 

< s irk, < A19 > 

and 

The tangential part at r = R reads for the inner solution 

Tj B = E (2(W^ 2 )(^ - I)*?™ + «r&) V s tt >m (*> 0), (A21) 
and for the outer solution 

Tf' = E (-2(V/R 2 )(l + 2)$S + arVS) V 5 3M#, 0), (A22) 

where 

vS "fii? e ^ + S eA (A23) 

Accomodation of these well-known properties of Lamb's solution to our problem requires 
to balance the stress discontinuities at the vesicle surface with those exerted by the vesicle's 
bending moments and its inhomogeneous surface tension. The key approximation involved 
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for quasi-spherical vesicles is to assume that all stresses act on a sphere of radius R rather 
than on the time-dependent vesicle's shape. This approximation can be controlled for small 
external flow and small excess area. However, in the absence of a workable alternative, we 
will use it for all flow strengths. The stress balance thus reads 

1 SF 

T out = T m + / A24 ) 

R R v^ R 



The stress exerted by the vesicle follows from (|TT|) in an expansion in spherical harmonics 
(k/R 2 ) ]T (EiU hm + 2a l , m )y lim {9, 0)e r - (T hm V s y l>m {0, <p)) . (A25) 



as 

1 5F 

Evaluating the tangential balance leads to 



21 ' rl -(X l>m + 2Y ljm ). (A26) 



1(1 + 1) 



Evaluating the normal balance, using this result for ai >m , and setting X Z m = Xf^ and 
Yl,m — Yfm leads to the relations (|37l) -(|39l) quoted in the main part. 
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